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ABSTRACT 

We analyze the time dependence of fluid variables in general relativistic, magnetohydrodynamic 
simulations of accretion flows onto a black hole with dimensionless spin parameter a/M = 0.9. We 
consider both the case where the angular momentum of the accretion material is aligned with the 
black hole spin axis (an untilted flow) and where it is misaligned by 15° (a tilted flow). In comparison 
to the untilted simulation, the tilted simulation exhibits a clear excess of inertial variability, that 
is, variability at frequencies below the local radial epicyclic frequency. We further study the radial 
structure of this inertial-like power by focusing on a radially extended band at 118(M/lOM0)~^Hz 
found in each of the three analyzed fluid variables. The three dimensional density structure at this 
frequency suggests that the power is a composite oscillation whose dominant components are an over 
dense clump corotating with the background flow, a low order inertial wave, and a low order inertial- 
acoustic wave. Our results provide preliminary confirmation of earlier suggestions that disk tilt can 
be an important excitation mechanism for inertial waves. 

Subject headings: accretion, accretion disks — black hole physics — MHD — turbulence — waves — 
X-rays: binaries 



; 1. INTRODUCTION 

Quasi-periodic oscillations (QPOs) are observed in the 
] X-r ay light curves of many black hole X-ray binaries 
(see [Remillard fc McClintockl l2QQ6l for a recent review) . 
They have also been o bserved in extragalactic u l tralu- 
' minous X-ray sources (|Strohmaver fc Mushotzkvl l2QQ3l : 
. [Strohmaver et al. 2007) and in one active galactic nu- 
cleus (Gierlihski et al. 2008). The origin of these phe- 
nomena is still far from clear. One class of models cen- 
ters on trapped wave modes within the accretion flow. 
In particular, a variety of modes have been proposed 
in hydrodynamic models of geometrically thin accretion 
disks (e.g. Wagoner 1999; Kato 2001). Axisymmetric 
inertial-acoustic modes ("/-" or "inner p-modes") and 
axisymmetric inertial modes ("^-" or "r-modes") ^ can 
be trapped due to the existence of the innermost sta- 
ble circular orbit around a black hole, which produces 
a maximum in the radial profile of the radial epicyclic 
frequency. Non- axisymmetric inertial modes can also be 
radially trapped between inner and outer Lindblad reso- 
nances. 

The very existence of inertial modes in accretion 
disks has been challenged recently, as numerical simu- 
lations appear to show that magnetorotational (MRI) 
turbulence suppresses them. This has been demon- 

^ Throughout this paper we will use the term inertial mode and 
r-mode interchangeably. We prefer to avoid the use of the term 
5f-mode, whose primary restoring force is generally due to entropy 
gradients. In disks, the primary restoring force for inertial modes 
is due to specific angular momentum gradients. 



strat ed in local shearing box simulations (I Arras et al.l 
120061), whose boundary conditions would artificially trap 
a discrete axisymmetric mode spectrum. Power spec- 
tra from these simulations show discrete acoustic modes 
and a radial epicyclic oscillation, but no inertial modes. 
This behavior has also been seen in global simula- 
tions of accretion disks in a pseudo-Newtonian poten- 
tial: hydrodynamic disks exhibit trapped axisymmetric 
inertial modes, whereas MHD turbulent disks do not 
(Reynolds fc Miller 2008). It appears that axisymmet- 
ric inertial modes, which necessarily have frequencies at 
or below the local radial epicyclic frequency, are partic- 
ularly vulnerable to nonlinear damping by MRI turbu- 
lence, which has a power spectrum that peaks near the 
orbital frequency. Non-axisymmetric inertial modes can 
have higher frequencies and might be less vulnerable to 
MRI turbulence, but there has as yet been no convinc- 
ing demonstration of the existence of a trapped inertial 
mode spectrum in any MHD simulation. 

Another reason why discrete inertial modes may not be 
present in accretion disks with MRI turbulence is that 
even subthermal magnetic fields can extend the inner 
trapping radius of the propagation zone of both axisym- 
metric and non-axisymmetric inertial m odes down to th e 
innermost stable circular orbit (ISCO) (|Fu fc Laill2009h . 
Unless inertial waves can reflect off the plunging region 
of the flow, standing waves will no longer be sustainable. 

On the other hand, it has been suggested that warps 
and eccentricity in disks may play a fundamental role 
in a nonlinear excitation mechanism of inertial modes 
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in geometrically thin ac cretion disks (|KatQ|[2QQ4al l2QQ8l : 
iFerreira fc Ogilvi3l2QQ8l ) . These large scale deformations 
interact with trapped r- modes giving rise to intermediate 
modes. These intermediate modes then couple back to 
the warp or eccentricity to produce positive feedba ck on 
the original r-mode oscillations (see Ferreira fc Qgilvid 
^008 for a more detailed explanation of this coupling 
mechanism) . 

[Fragile et alJ (|2QQ7l ) recently completed fully general 
relativistic, MHD simulations of accretion disks with 
misaligned black hole spin and disk angular momen- 
tum vectors ("tilted disks"). These tilted disks are 
strongly warped near the black hole and exhibit global 
epicyclic oscillations superimp osed on the turbulence in 
the flow (Fragile fc Blaesl[2QQ8[ ). These oscillations mani- 
fest themselves as eccentric orbits of fluid elements in the 
disk, albeit with a 180° flip in orientation of the elliptical 
orbits across the midplane of the disk. Tilted disks may 
therefore provide favorable conditions for the nonlinear 
excitation of inertial modes. 

Here we report the results of a search for trapped 
modes in simulations of both an untilted and tilted disk. 
In agreement with previous work, we do not observe the 
presence of modes in the untilted simulation. However, 
at the same numerical resolution, we observe signiflcant 
excess power with frequencies characteristic of inertial 
waves in the tilted disk. Because of the complexity of 
the background flow, the physical nature of this power is 
difficult to determine. Nevertheless, its spatial structure 
appears to conflrm that this is partly inertial in charac- 
ter. This may represent preliminary conflrmation of the 
warp/eccentricity excitation mechanism, showing that it 
can be strong enough to overcome damping due to MRI 
turbulence. It also shows that trapping, in the sense of 
radial localization of power, can be maintained even in 
the presence of magnetic flelds. 

This paper is organized in the following manner. In 
section 2 we briefly review the simulation parameters and 
discuss our power spectrum analysis of the simulation 
data. Section 3 presents power spectra from both un- 
tilted and tilted geometries, evidence that the latter may 
be dominated by inertial-like variability, and an analy- 
sis of the three dimensional structure of the power at 
a particular frequency within this inertial regime. We 
summarize our conclusions in section 4. 

2. SIMULATIONS AND TIME SERIES ANALYSIS 

We analyz e the two hig;h resolution simul ations of 
[Fragile et all (|2QQ7D and [Fragile fc BlaesI (|2QQ8D : an un- 
tilted simulation (90h) with the initial fluid angular mo- 
mentum aligned with the spin axis of the Kerr black hole, 
and a simulation (915h) with the initial fluid angular mo- 
mentum tilted by 15° relative to the spin axis of the hole. 
The black hole spacetime in each simulation is described 
using Kerr-Schild coordinates, and has a spin parame- 
ter a/M = 0.9. Apart from the tilt, both simulations 
start with identical initial conditions: a torus with pres- 
sure maximum centered at 25i?G, where Rq = GMjc? ^ 
seeded with weak poloidal magnetic fleld loops that fol- 
low the equipressure surfaces. Both simulations evolve 
the disk through ten test particle orbital periods calcu- 
lated at the initial pressure maximum (a time unit here- 
after generically referred to as an "orbit"). 

Again apart from the tilt, these simulations run on es- 



sentially the same nested grids with equivalent peak res- 
olution of 128^ zones realized throughout the bulk of the 
simulation domain. These grids have suflicient spatial 
resolution at all radii t o reasonably capture the MRI. As 
discussed in Stone et al. (1996), several grid zones are re- 
quired per critical MRI wavelength (Amri = 27rvA/^ = 

^ I ^.ix pVfi ^ where is the Alfven velocity, is 
the square of the magnetic fleld, p is the fluid mass den- 
sity, and Vt is the local orbital frequency) to ensure a 
convergent treatment of the instability and associated 
turbulence. As illustrated in flgurelH there are at least 
10 grid zones per critical wavelength, apart from small, 
isolated, transient patches where the fleld is weak. 

We restrict our analysis to a poloidal sector of grid 
zones within 35° of the original outer torus midplane at 
all radii and azimuthal angles. Within this region, we 
uniformly sample fluid variables at each grid zone 200 
times per orbit (that is, 7 times per orbital period of the 
direct ISCO on the black hole equatorial plane). Our 
timing analysis begins four orbital periods after the start 
of each simulation to ensure fully developed accretion and 
MRI turbulence, then extends for six additional orbits. 

We study the time variation of three fluid variables, 
namely, the rest mass density p as measured in the fluid's 
rest frame and the radial and poloidal components of the 
3-velocity: = /u^ and = /u^ . Here, are the 
general relativistic 4-velocity components expressed in 
tilted Kerr-Schild coordinates (Fragile et al. 2007), mea- 
sured at zone centers. 

We search for quasi-periodic activity in these vari- 
ables through analysis of radially dependent temporal 
power spectra. The three variables are measured in each 
grid zone at regular time intervals, and linear secular 
changes are subtracted from the zone's time series as 
in Schnitt man et al.i (;2006i) . Other authors have used 
more elaborate and physically-motivated fltting func- 
tions to accomplish this "pre- whitening." For example, 
iRevnolds fc Millerl (|2008l ) advocate exponential forms 
because of the secular spreading of mass within the disk 
which, in our simulation and depending on the location 
of a given cell with respect to the initial torus, can act 
to increase or decrease density. On the other hand, in 
a tilted conflguration, another signiflcant secular change 
is Lense-Thirring precession with a period of roughly 90 
orbits. In the six analyzed orbits, the disk's angular mo- 
mentum vector precesses roughly 24° about the black 
hole's spin axis. This precession induces some fraction 
of the full period of a sinusoidal modulation at the pre- 
cession frequency. Whether this results in an increase or 
decrease in density is determined by the location of the 
given zone with respect to the initial torus. Since these 
two effects are completely independent of each other and 
can act to either increase or decrease the local density, it 
is possible for the two to act in opposition. We employ, 
therefore, a simpler linear (i.e. model-independent) fit 
to remove secular changes. Finally, these pre- whitened 
time series are windowed using a Bartlett window, zero 
padded and fast Fourier transforni ed to produce the 
zones' raw spectra (jPress et al.lll992l ). 

In calculating the final radially dependent, temporal 
power spectrum, we average our spectra on shells of con- 
stant coordinate radius and divide by appropriate back- 
ground scaling factors to enhance the radial contrast be- 
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Fig. 1. — Number of azimuthal grid zones per critical MRI wavelength in the torus' midplane for the untilted (left) and tilted (right) 
simulations. Excepting a few small, isolated, transient patches, the entire simulation domain has roughly 10 or more zones per wavelength, 
indicating MRI turbulence should be well captured. 



tween features using physically relevant scales. Symbol- 
ically, 

Pp{rJ)^{\p\')e,^/{p')tA^ 
P.ArJ)^{\vr\')e,^/{cs)le,^ (1) 

Pvo f) = (l^~d^)^,0/(Cs)t%,0, 

where p, Vr^ and vq are the pre- whitened, windowed, 
padded, and Fourier transformed fluid variables as func- 
tions of space and frequency. The weight factors {p^)t,e,(t) 
and {cs)t^Q^^ are the mean square density and squared 

average sound speed, Cg = y^P/p^ respectively. The ra- 
dial profiles of these weight factors are shown in figure 
[21 We choose the sound speed as an appropriate scaling 
speed for both acoustic and inertial modes because both 
mode types involve fluid oscillations with vertical struc- 
ture of order the disk scale height, and frequencies 
comparable to the orbital frequency, Q. Thus, the rele- 
vant velocity scale is HQ which, by vertical hydrostatic 
equilibrium, is approximately equal to the sound speed. 

3. RESULTS 

Figure [3] depicts the radial power spectrum in the un- 
tilted and tilted simulations. In the untilted simulation, 
three spectral features are prominent. First, high power 
at large radii in all variables reflects artificially enhanced 
variability by the extremely low background density and 
sound speed weight factors in those regions (see figure [2|). 
Very little material is present in these regions since the 
simulated time period is short when compared with the 
diffusion time required to spread the disk to such large 
radii. 

Second, density variability at high frequencies near and 
within the IS CO is present. We do not fully understand 
this behavior. Similar variability seems absent in this 
simulation's velocity variables and, as will be discussed 



later, in all three variables in the tilted simulation. This 
suggests that the feature is physical (i.e. it is not a relic 
of the Fourier transform procedure) and reflects a funda- 
mental difference in the accretion flow between the two 
disk configurations. For a detailed comparison of the ti- 
tled and untilted flow patterns, see Fragile et al. (2003)- 
Third, significant power tracks the orbital frequency 
(or, arguably, th e vertical epicyclic f r equen cy) in the den- 
sity spectrum. iRevnolds fc Millerl (|2QQ8l ) suggest that 
this power is a result of local vertical epicyclic oscilla- 
tions and would predict higher order acoustic oscillations 
with frequencies 



uj :^ [n(n7 — n 3 — 7)/2]"^/^/^6i, 



(2) 



where n = 1,2,3,... is a vertical quantum number, 7 is 
the adiabatic index (which is 5/3 in o ur simulations), an d 
K0 is the vertical epicyclic frequency (jBlaes et al.ll2QQ6l ).^ 
While such v ertical acoustic mod e s aptly explain the 
features seen in IRevnolds fc Milled ' s (|2QQ8l ) simulations, 
they do not appear to correctly describe our spectra. 
Even though we cannot claim the power more closely 
tracks the orbital frequency than the vertical epicyclic 
frequency, it appears that the higher order diagonal 
tracks tightly follow harmonics of the orbital frequency 
with integer coefficients (indicated by dotted lines in fig- 
ure [3]) rather than those curves described by equation [2] 
relating to the vertical epicyclic frequency. In particular, 
we note this agreement even at frequencies higher than 
the vertical epicyclic frequency maximum. Additionally, 

^ Note that k,q ^ fl for our Kerr spacetime, in contrast to th e 
pseudo-Newtonian spacetime used by Revnolds Mille r^ (I 2008I V 
Their simulations did, in fact, exhibit modes with frequencies 
u; ~ (727 + ly^^hiQ with n = 0,1,2,..., because they employed 
an isothermal vertical profile with an adiabatic equation of state. 
Our simulations are better described by an adiabatic profile and 
equation of state, thus acoustic waves should obey the dispersion 
relation [2] 
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Fig. 2. — Normalization curves as a function of coordinate radius r showing the mean squared density (top) and squared mean sound 
speed (bottom) for both the untilted (left) and tilted (right) simulations, in simulation units (i.e. the natural density po = /G^M'^ and 
the speed of light c). 




Fig. 3. — Shell averaged power spectral density as a function of frequency / and coordinate radius r in density (top), radial velocity 
component Vr (middle), and polar velocity component vq (bottom), for the untilted (left) and tilted (right) simulations. Overplots include 
the orbital frequency (solid) and its harmonics (dotted), the geodesic radial epicyclic frequency and the ISCO radius (dashed), and the 
geodesic vertical epicyclic frequency (triple-dot dashed). Note that, for simplicity, an arbitrary constant has been subtracted from each 
logarithmic plot, setting the colorbars' centers at 1.5. 
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tracks do not appear strongly in the polar velocity spec- 
trum as would be expected if acoustic motion in that 
direction were present. Hence, we argue that the funda- 
mental dynamics giving rise to our tracks is not related 
to vertical acoustic motion. 

An alternative interpretation for the power along the 
orbital frequency and its harmonics in our simulations 
is that they represent non-axisymmetric clumps orbit- 
ing the central black hole on Keplerian trajectories. One 
would not expect the presence of such orbiting clumps 
in th e [Reynolds & Miller (2008) simulations (Reynolds 
[2009). They employed a 30° azimuthal wedge with pe- 
riodic boundary conditions in ^, so their lowest order, 
non-axisymmetric structure would, in our language, have 
at least an azimuthal quantum number m = 12. These 
high m modes would manifest as fluid oscillations at a fre- 
quency 12Q but woul d likel y be dominated by turbulence 
in both Revnol ds fc Miller! ' s f2008) simulations and our 
own. In addition, it could be that the low m clumps that 
we see in our simulations dominate the vertical acoustic 
waves that they saw in their wedge. 

The tilted simulation's power spectrum, shown in fig- 
ure [3l similarly exhibits power along the orbital fre- 
quency and its harmonics and spurious variability at 
large radii. In contrast to the untilted configuration's 
spectra, no high frequency variability appears near or 
within the IS CO in any variable. A characteristically dif- 
ferent shape to the overall spectra is evident: while the 
untilted spectra are dominated by acoustic-like variabil- 
ity (i.e. variability at frequencies larger than the orbital 
frequency), the tilted spectra exhibit considerably more 
power at inertial-like frequencies (i.e. frequencies smaller 
than the orbital frequency, particularly those bounded 
above by the radial epicyclic frequency) . To visualize the 
stark contrast between the the untilted and tilted power 
spectra, figure [4] displays the power integrated in radius 
between i^isco (~ 2.4i^G) and lOi^c- It is clear that the 
tilted simulation is contrastingly dominated by power at 
frequencies characteristic of inertial variability, whereas 
the untilted simulation is not. We believe this power has 
not been previously observed in MHD simulations. 

Figures [3] and [4] may also indicate higher power peaks 
amid the continuum of inertial-like variability. Localized 
in frequency and extended in radius, such peaks suggest 
coherent variability arguably spanning the region from 
the spacetime's ISCO to the radius at which the orbital 
frequency equals the oscillation frequency. Additionally, 
power seems to leak inwards of the ISCO, particularly in 
the density spectrum. We can identify matching peaks in 
each of the three fluid variables at corresponding frequen- 
cies, suggesting that the features are physical in nature 
and are not numerical artifacts. Among these, power 
at 118Hz is elevated above the continuum in the den- 
sity and radial velocity spectra and exhibits a minimal 
though identifiable excess in the polar velocity spectrum. 
We hope to characterize the nature of the general power 
at inertial frequencies by carefully analyzing the three 
dimensional structure of the variability at 118Hz. 

Isothermal diskoseismic models of untilted, flat disks 
indicate that perturbations with l ocal radial wave num - 
ber k obey the dispersion relation (|Okazaki et al.|[T987l ). 

2 _ {u;^-f^l){u^^-nf^l) 



where cu = co — mVt is the Doppler-shifted wave fre- 
quency and Kjf IS the radial epicyclic frequency.^ Wave 
propagation is only possible where /c^ > and equa- 
tion [3] is quadratic in o)^, so it is clear that two dis- 
tinct oscillation modes exist for a given frequency uo 
and set of vertical and azimuthal wave numbers n = 
0, 1, 2 . . . and m = 0, 1, 2 . . ., re spectively. As summa- 
rized in lFerreira fc Ogilvid (|2008l ). modes that have pres- 
sure as their dominant restoring force (p-modes) can 
exist at radii where uj < mQ — max(A>:r, V^^^*) and 
CO > mr^ + max(/^^, ^/n^^e)•, while modes whose dominant 
restoring force is inertial in character (r-modes) can lie 
where mVt — min(/^^, \/nf<io) < uo < mVt + min(/^^, ^/nno). 
Outside a given mode's trapping region (i.e. where k 
is imaginary), the oscillation damps exponentially. Un- 
fortunately, this diskoseismic theory applies to perturba- 
tions of a geometrically thin, flat accretion disk. It is 
therefore difficult to compare its predictions directly to 
a pressure supported, moderately tilted disk due to the 
nonlinear effects of spatially varying tilt and twist and 
significant vertical thickness. 

Power at 118Hz displays complex radial structure. 
Two local maxima appear just inside the inner and outer 
edges of the radial epicyclic frequency, and another lies 
just within the orbital frequency's corotation radius. The 
first two features may imply a superposition of modes 
trapped within respective resonances, while the third 
has the characteristics of the aforementioned Keplerian 
clumps. There is also significant mode power well inside 
the equatorial test particle ISCO, which might be due 
to several factors. Pressure gradients can sustain sta- 
ble fluid element orbits down to a smaller effective ISCO 
radius. Additionally, mode power may manifest itself in- 
side this effective ISCO due to active accretion; any mode 
with a trapping radius close to the plunging region may 
regularly slough off material to smaller radii at the mode 
frequency. Such fluid pulses would reveal themselves in 
spectra as higher power at all radii within the plunging 
region. 

To quantify the radial structure and better understand 
the physical nature of the 118Hz variability, we extract 
its spatial structure from the real and imaginary parts of 
the Fourier transform of our fluid variables. Diskoseismic 
eigenmodes possess well defined azimuthal wave numbers 
m and parity with respect to the disk equatorial plane. 
Because a tilted disk lacks both axisymmetry and verti- 
cal reflection symmetry, its wave eigenmodes need not be 
characterized by such quantum numbers. Still, in order 
to connect to the analytic theory of untilted disks, we 
project the 118Hz density structure onto six orthogonal 
basis functions with well-defined m and parity. Figure [5] 
depicts 3D isosurfaces of the real part of the eigenfunc- 
tion and its projections; an animation of these isosurfaces 
may be found online with this publication. A clear m = 1 
structure with even parity is the dominant component, 
accounting for roughly 50% of the total power. The next 

3 lOkazaki et"aD ([1983) assume an isothermal equation of state 
and vertical profile. For an adiabatic equation of state and vertical 
profile, the acoustic factor must be slightly modified as in equation 
[2] above. However, for the important n = and n = 1 modes, both 
dispersion relations agree. Physically, these two modes exhibit no 
vertical oscillation and a bulk vertical oscillation of the local center 
of mass, respectively. There is therefore no internal compression, 
and the thermodynamics are irrelevant. 
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Fig. 4. — Untilted (left) and tilted (right) disk's power spectra averaged radially between i^isco and IORq and on shells as a function 
of frequency / in density p (top), radial velocity Vr (middle), and polar velocity vq (bottom). The curves' fundamentally different shapes 
indicate that while the untilted simulation is dominated by acoustic-like variability, the tilted simulation is better characterized by inertial- 
like spectral power. 



two most powerful components are the m = and m = 2 
modes with odd parity, possessing roughly 10% and 7% 
of the total power, respectively. 

The radial distribution of power in each of these pro- 
jected modes sheds light on the physical nature of the 
oscillations. Figure [6] illustrates the total 118Hz power 
and each of its three most dominant mode components as 
a function of radius. Additionally, vertical dotted lines 
indicate the zeros of /c^ for each mode's quantum num- 
bers and horizontal lines show the regions in which k'^ 
is positive, that is, the mode propagation regions. For 
the (odd, m = 1) and (odd, m = 2) modes, local max- 
ima occur within predicted trapping regions, suggesting 
that where the 118Hz feature is most significant, the for- 
mer behaves as an r-mode bound by the radial epicyclic 
frequency, and the latter is a p-mode trapped between 
2Q — tvQ and a reflecting boundary near the plunging re- 
gion. Conversely, the (even, m = 1) mode peaks near the 
orbital frequency (indicated by the vertical, black dotted 
line) and well inside a region forbidden by the disko- 
seismic models. Recalling the trail of high power peaks 
tracking the orbital frequency in figure [3l we speculate 
that this (even, m = 1) component is an orbiting clump 
and not a trapped inertial or acoustic mode. 

However, we note an alternative explanation for power 
in the (odd, m = 0) and (odd, m = 2) components. 
Consider the following: suppose we have a pure (even, 
m = 1) mode but that its rotation axis is slightly mis- 
aligned with the z-axis of a cylindrical coordinate sys- 



tem. If we decompose this configuration with respect to 
the true z-axis, we will find some small amount of power 
has spilled into the other modes, most notably into the 
(odd, m = 0) and (odd, m = 2) modes. In a prac- 
tice, it is difficult to accurately determine the position of 
the disk midplane, and it is possible our decomposition 
mixes power between modes. However, the correlations 
between the predicted mode trapping regions and the ra- 
dial power profile suggest the power distribution is not 
entirely a numerical artifact. 

In addition to the nonlinear excitation of trapped 
waves by warps and eccentric orbits within the disk flow, 
a titled accretion flow may be more conducive to mode 
trapping because of the distinct character of its plung- 
ing region. In contrast with an untilted disk which has 
a largely axisymmetric plunging region, mass accretion 
at small radii in our tilted simulation is confined to two 
discrete plunging streams f Fragile e t al.ll2QQ7i ). Thus, it 
is conceivable that the azimuthal sections of the inner 
disk away from these streams provide significant radial 
reflection of traveling waves in the disk. 

4. CONCLUSIONS 

Previous work has demonstrated that inertial modes 
are vulnerable to disruption by MR I turbulence 
(| Arras et al.ll2QQ6l : lRevnolds fc Millerll2QQ8[ ). Our Fourier 
analysis of the untilted simulation, a configuration devoid 
of tilt or eccentricity, corroborates this conclusion. How- 
ever, nonlinear analytic calculations (^Kato..2QQ4al l2QQ8l : 
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Fig. 5. — Positive (red) and negative (blue) isosurfaces of density fluctuations evaluated at 118Hz in the tilted disk simulation. Subflgures 
include the unmodifled complete mode structure (leftmost top), projections onto six orthogonal basis functions having even or odd parity in 
the polar direction and azimuthal quantum numbers m = 0, 1, or 2 (right three columns), and the residual, that is, the difference between 
the complete mode and the sum of the six projections (leftmost bottom). The gray axes represent a rotated coordinate system in which the 
z-axis is normal to the average orientation of the disk midplane. Parity is deflned with respect to this xy-plane. A green vector indicates 
the black hole spin axis. Also given is the fraction of the total mode power in each projection. The surface levels are drawn at half the 
maximum of the unmodifled mode's absolute value times the square root of the respective projection's fractional power (or simply times 
1 in the unmodifled case). While this static plot depicts the real (or zero phase) component of the eigenmode, the online animated flgure 
varies the mode's complex phase as a function of time. 



iFerreira fc Ogilvi3l2QQ8l ) have indicated that disk warps 
and eccentricity may excite inertial modes; both of these 
features are generically present in tilted accretion flows. 
Though the analysis of the variability in our tilted sim- 
ulation at 118Hz does not appear to be wholly inertial 
(or even acoustic) in nature, the presence of two weak, 
odd parity, m = and m = 2 components provides some 
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Fig. 6. — Radial distribution of power at 118Hz (black), its 
(even, m = 1) component (green), (odd, m = 0) component (red), 
(odd, m = 2) component (blue), and residual (dashed black). Col- 
ored, vertical, and dotted lines indicate the zeros of the squared 
radial wave number for each mode, where even parity has been 
naively associated with n = and odd with n = 1. Horizontal lines 
designate the regions in which is positive, that is, the regions 
permitting wave propagation according to equation [3] The orbital 
frequency is indicated by a black, vertical, dotted line. 



support for this hypothesis. 

Evidence that disk tilt is connected with the excitation 
of inertial or acoustic waves is twofold. First, fundamen- 
tally different shapes in the spectra from the untilted and 
tilted simulations implies a correlation between tilt and 
the power at frequencies characteristic of inertial waves. 
Second, the radial structure of the variability suggests 
a superposition of modes, at least one of which may be 
described as inertial in nature in the context of relevant 
analytic models. Together, these clues establish the pres- 
ence of variability unique to a tilted geometry which ap- 
pears to be at least partially inertial in character. 

Recalling that a tilted disk configuration naturally 
yields a low frequency QPO activity in black hole X-ray 
binaries and speculating that our work's inertial-like vari- 
ability may be related to high frequency QPO activity in 
the same systems, we advocate further numerical explo- 
ration of these tilted geometries. Torques associated with 
the Kerr spacetime can produce a bodily precession of the 
tilted disk at frequencies within the 0.1 and 30Hz range 
associated with the low frequency QPO. Because this 
precession frequency is strongly dependent on the radial 
extent of the disk in this model, it readily explains the ob- 
served correlation between the LFQPO's frequency and 
disk fiux tRemihard^fe McClintoc k 2006; Fragile fc Blaesj 
120081 : llngram et al. I I2009D . Shnilarlv. our current work 
shows inertial or acoustic variability excited and main- 
tained by the disk tilt that exhibits frequencies of order 
the radial epicyclic frequency maximum. Such frequen- 
cies fall within the observed range of 40 to 450Hz for 
HFQPOs in stellar mass black hole systems. 
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HFQPOs are only observed when a system's X-ray 
emission exhibi ts the characteristics of the ste ep power- 
law state (cf. iRemihard fc McClintocg I2QQ6I for a re- 
view of the properties of X-ray black hole binaries), 
and we must acknowledge, however, that our simula- 
tions are almost certainly not accurate representations 
of that state. For instance, since they neither account 
for radiative cooling nor capture all forms of dissipa- 
tion, our simulations likely do not reproduce some of 
the key features found in various models for the state's 
structure (e.g. the energetically significant corona in 
[Done, Gi erlihski. fc Kubota I [2007). However, we are 
confident that our simulations and analysis provide in- 
sight into the dynamical effects of a tilted accretion 
fiow, particularly, the excitation of inertial and acoustic 



modes. 
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